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Abstract 

There is proposed a model of scale-free random graphs which are 
locally close to the uncorrelated complex random networks with diver- 
gent (fc 2 ) studied in e.g. S. N. Dorogovtsev et al, Rev. Mod. Phys. 
80, 1275 (2008). It is shown that the Ising model on these graphs 
with interaction intensities of arbitrary signs with probability one is in 
a paramagnetic state at sufficiently high finite values of the tempera- 
ture. For the same graphs, the bond percolation model with probability 
one is in a nonpercolative state for positive values of the percolation 
probability. Possible extensions are discussed. 

1 Introduction 

In statistical physics, a paradigm model for cooperative phenomena is the 
Ising model of interacting spins attached to the vertices of a graph G = 
(V, E). This model is described by the Hamiltonian 

H = — JxyVxPy — h'^dx, (1.1) 

{x,y}eE xev 

where J xy are interaction intensities, h is an external field, and the spins 
a x take values ±1. The sums run over the sets of edges E and vertices V 
of the graph G, which is supposed to be connected and countably infinite. 
In view of the latter fact, the Hamiltonian Ql.ip has no direct mathematical 
meaning and is used as a formula for local Hamiltonians, H\, which one 
obtains by restricting the sums to finite subsets A C V and E/\ C E, where 
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Ea is the collection of edges with both ends in A. Then the thermodynamic 
properties of (jl.ip are studied in the limit A — > V. In a more sophisticated 
theory [I], a thermodynamic state is a Gibbs measure, and a phase transition 
is associated with the possibility for multiple thermodynamic states to exist 
at the same values of the temperature and h. 

Obviously, the properties of the model (jl.ip are closely related to those 
of the underlying graph. The case of a special interest is where this graph 
is random, with vertices having arbitrary number of neighbors. Then the 
graph itself can be characterized by the properties of the corresponding 
Ising model. Perhaps in view of this fact, methods of statistical physics are 
widely used also in the study of the underlying random graphs. A long list 
of publications on this topic can be found in the review articles [2j [31 B] • 

For a vertex x, the degree k(x) is defined as the number of neighbors of x, 
i.e. the number of y £ V such that {x, y} € E. If G is random, the degree of 
each x € V is an integer valued random variable and hence is characterized 
by the probability distribution 

Prob{/c(x) = k} = p k (x), k > 1, x e V. (1.2) 

In simple cases, the degrees k(x) are independent and identically distributed, 
which in particular means that the probabilities Pk(x) are the same for all 
x. If pk = Ck~ x , at least for big enough k and an appropriate C > 0, the 
corresponding graph is said to be scale- free, c.f. [5]. In this case, the m-th 
moment 

{k m )=Y,k m p k , m>l, 

exists provided A > m + 1. Notably, (A; 2 ) = +oo for many real- world net- 
works [1]. 

An example of a random graph with independent k(x) is a Galton- 
Watson tree, which is the genealogic tree of a Galton- Watson branching 
process [SI Ej. In more complex models, the graph G is obtained, more or 
less explicitly, as the limit of an increasing sequence of finite random graphs, 
(V n , E„), n 6 N, see, e.g. [HJ [TOj, [11] . In the physical literature, such con- 
structions are rather informal that leaves certain possibilities for different 
interpretations, which we are going to use in this work. In 2002, almost 
simultaneously there appeared two papers [121 H3] where the model (jl.ip on 
a scale-free graph was studied. The starting point is a scale-free distribution 
{PfcjfceN, which "characterizes the entire net", whereas the "local structure 
of connections" is characterized by "the distribution of the number of con- 
nections of the nearest neighbor of a vertex", which is taken to be kpk/{k), 
see the beginning of the first section in [12]. In [13], a similar assumption 
was made, see equation (3) thereof, and also (1) in [3]. In the latter article, 
such graphs are called complex networks in contrast to the classical models 
called simple networks, one of which is presumably the Galton- Watson tree 
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mentioned above. In both papers [X2^ 113]. it was found that the model (jl.ip 
with h = 0, ferromagnetic interactions J xy = J > 0, and (k 2 ) = +oo is in 
an ordered (ferromagnetic) state at all temperatures^ T > 0. Recently, in 
[10\ [TT] the mathematical construction of a random graph model was pre- 
sented, in which the mentioned properties were realized. For this model, the 
conclusion of \12\ \13\ can be mathematically proven, see [15] and Subsection 
15.21 below. The aim of this work is to construct "uncorrelated complex net- 
works" with (k 2 ) = oo, consistent with the assumptions made in [3l 112} [T3], 
for which the model (jl.ip with interaction intensities of arbitrary signs can 
be in a paramagnetic state at finite temperatures. In Section [21 we perform 
this construction and obtain a graph model, which is a combination of the 
one used in [10l E] and of the ordinary Galton- Watson tree. In Section [3l 
we prove Theorem 13.11 the main result of which is that the model (jl.ip on 
the introduced graphs, with any signs of J xy , can be in a paramagnetic state 
at finite temperatures and h = 0, provided (k 1+a ) < oo for some a > 0. 
This covers the case of scale-free distributions with A G (2,3). The proof 
of Theorem 13.11 relies upon two lemmas, proven in Section SJ The result of 
Theorem 13.11 and the graph model introduced in Section [2] are discussed in 
the concluding Section El where we give also an extension of the results of 
[12\ [P3] concerning the model (|1.1|) with interaction intensities of arbitrary 
signs. Except possibly for Section [H the paper is reasonably self-contained 
and accessible to non- mathematicians. 



2 The graph models 

Since the Galton- Watson tree will serve us as a prototype, we beging by 
presenting this model. More on this topic can be found in \7\ \W\. 

2.1 The Galton- Watson tree 

To exclude the possibility of obtaining a finite tree, or a tree of bounded 
degree, we shall suppose that the basic random variable X takes positive 
integer values k G N with probability Vroh{X = k) = pt+i such that 

n 

p k G [0, 1), J^pk = 1, < 1, for any n G N, (2.1) 

k>2 k=2 

a := kPk < °°- 

k>2 

One observes that 

a > 2, (2.2) 

1 In contrast to regular graphs, e.g. crystal lattices, the Ising model on sufficiently 
'dense' non-random graphs can be in an ordered state at all temperatures, see the example 
in [H]. 
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which readily follows form the properties of p assumed in (|2.ip . Set Lq = 1, 
L\ = X, and 

L n -i 

L n = 'y ] X n j, (2-3) 
i=l 

where all X n j are independent copies of X, see Chapters II and III in [6] or 
Chapter 5 in [7 J . Then L n is the number of individuals in generation n of a 
Galton- Watson branching process. The process starts with one individual, 
and each individual in generation n produces offsprings at random, inde- 
pendently and with the same distribution as X. The individuals existing 
in all generations form the vertex set V of the corresponding (genealogic) 
Galton- Watson graph G, with undirected edges {x,y} € E connecting each 
individual and its offsprings. Under the assumptions made in (|2.ip . the cor- 
responding process does not extinct and L n —¥ +oo with probability one. 
Hence, G is an infinite tree - an acyclic connected graph, which has no 
leaves, except possibly for the root o, the degree of which coincides with the 
number of the corresponding offsprings Lq. The degree of any other vertex 
equals the number of its offsprings plus one. By construction, the vertex 
degrees are independent identically distributed random variables such that 

Prob{fc(o) = k} = Pk+i, Prob{k(x) = k} = p k , xGV\{o}. (2.4) 

Now let us fix the graph-theoretical and probabilistic terminology used in 
this article. For a graph G = (V, E) and vertices x, y E V, a path y) 
is a sequence {lo^i, • • • ,x n }, n € N, such that: (a) xq = x, x n = y; (b) 
{xi, xi + i} € E; (c) neither of xq, . . . , x n -\ can be repeated. The length of the 
path is ||i?(x,y)|| = n; x is its origin and y is terminus. A ray with origin xo 
is an infinite sequence of distinct vertices {xq, x\, . . . , x n , . . .} with property 
(b). A graph ensemble Q is the collection of all possible realizations, e.g., 
as just described. A random graph model is a triple {Q,F,P), where T is 
a certain family of subsets of Q and P : J- — > [0, 1] is probability. For the 
Galton- Watson tree, P is defined by the collection p = {pk}- For A € J~, 
P(A) is the probability of the event: "the graph chosen at random from Q 
belongs to A". Some events can occur with probability either zero or one 
- the celebrated zero-one law. In the Galton- Watson tree, if A is the set of 
all realizations with k(o) = 1, then P(A) = p\. We say that a property of 
a random graph holds with probability one if this property is possessed by 
all of the elements of some A such that P{A) = 1. Note that this A need 
not be the whole ensemble Q. Usually, when one deals with models with 
independent vertex degrees, the triple (Q,F,P) does not appear explicitly. 
This will also be the case in the rest of this paper. 
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2.2 The configuration model 

Let p := {pk}k>2 be as in (|2.ip . and 

Pk ■= kpk/a, k>2. (2.5) 

"Uncorrelated complex networks" were 'defined' in [31 [121 E] by means of 
the following conditions. The main one is that "the neighbor of a vertex 
has degree k with probability pk" , which can be interpreted as the condition 
that in every path {xq, x\, . . . , x n }, the probability that k{x{) takes value k 
is pk- For convenience, we shall call it a path property. Another important 
assumptions are: (a) all k(x) are "uncorrelated" (in fact, independent); 
(b) two paths {x,xi, . . . ,y} and {x, x[, . . . , z}, with x\ / x[ and y ^ z, 
never intersect each other, i.e., the graph is a These assumptions 

are consistent with the following mathematical models. The first one is the 
ordinary Galton- Watson tree with p replaced by p, in [10] it is T(p, oo), see 
page 567 threof. A slightly different version is the so called configuration 
model, c.f. page 1279 in [3]. This model is a rooted tree with independent 
vertex degrees such that 

Prob(/c(o) = k) = pk, Pvob(k(x) = k) = pk- (2-6) 

This modef^ was studied in [TO] [11], in [10] it appears as T(P, p, oo). The 
main characteristic feature of both T(P, p, oo) and T(p, oo) is that each path 
x 7^ o, has the path property described above. The inconsistency 
with the root degrees does not change the global behavior of the graph, and 
hence of the corresponding model ([Lip . 

The construction of the configuration model can be visualized as follows. 
Consider a countable set of 'points' (a configuration), each of which is given 
a random mar 1(0 k £ N with probability distribution p as in (|2.ip . The marks 
of different points are independent. Then the graph is obtained as follows. 
The root degree is drawn from N according to Prob(fc(o) = k) = p^- If the 
result is m, the root receives m neighbors - 'particles' independently drawn 
from the configuration with probability proportional to the particle mark, 
i.e., Ck. Then the (unconditional) probability of drawing some particle, 
which bears mark k, is Ckp^ with C = 1/a since the total probability 
should sum up to one. In the result, each root neighbor has degree k with 
probability pk- In the same way one picks the next-neighbors of o keeping 
in mind that the only common element of the neighborhoods of different 
vertices is the root itself. This procedure is continued ad infinitum until the 
whole tree is constructed. An important observation here is that, except for 
the root degree, the distribution of k(x) is given by p, whereas p itself is the 
distribution of marks - auxiliary objects used in the graph construction. 

2 In [31 [T5] , the graphs were called "tree- like" if the intersections are negligible. Math- 
ematical definitions of tree-like graphs can be found in [101 1111 [15] . 

3 Also models which locally converge to T(P, p, oo), cf. Definition 2.1 in [10] . 
4 According to [3], a mark is the number of stubs. 
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2.3 The proposed model 

Here the main object of our study in this article is introduced. This is a 
random tree model, also consistent with the assumptions of O \12\ 113] . As 
in the case of the Galton- Watson tree, we beging with the description of the 
corresponding branching process. Let p, p, and a be as in i\2.1\\ and (|2.5|) . 
An additional parameter is s E N. The whole population falls into two 
types - 'distinguished' and 'ordinary' individuals. At the beginning, ordi- 
nary individuals are absent and there is one distinguished individual, which 
independently produces m offsprings with probability p m . If m < s, all of 
them are set to be distinguished. For m > s, s offsprings are set to be distin- 
guished and the remaining m — s ones are set to be ordinary. Each ordinary 
individual produces independently k ordinary individuals with probability 
Pk+i- Each distinguished individual, other than the initial one, produces 
independently m offsprings with probability p m +i- As above, if m < s, 
then all of the offsprings are distinguished. For m > s, s offsprings are 
distinguished and the remaining m — s ones are ordinary. This process is 
repeated ad infinitum. For n € No, by L n and L n we denote the number 
of distinguished and ordinary individuals in generation n, respectively. Let 
also X and X be independent and such that Prob(A = m) = p m +i and 
Prob(X = m) = p m +i, m 6 N. Then we have, cf. (|2.3|) . 



3=1 



(2.7) 



Ln—1 L n —\ 

i=i j=i 

where, for an integer x, we write x + = max{0, x}, and X n j and X n j are 
independent copies of X and X, respectively. From the point of view of 
ordinary particles, those represented by the second term in the expression 
for L n are immigrants, and their number in generation n is 

L n -1 

Y n =Y J (Xn,j ~ s) . (2.8) 
i=i + 

Then the total number of individuals in generation n can be written in the 
form 

Ln—1 L n —l 

L n := L n + L n = X n j + X n j. (2-9) 
3=1 3=1 

Since each distinguished individual can have at most s distinguished off- 
springs, we have that 

L n < s n . (2.10) 
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The graph in question is the genealogic tree of the process just described. 
Its construction can be visualized as follows. The graph is a rooted tree with 
Prob(fc(o) = k) = pk+i- Let k{6) = m < s. Then we draw m neighbors of 
o as in the configuration model. Denote them x\,... ,x m ; by construction 
Prob(/c(x;) = k) = pk for all I = 1, . . . m. If m > s, we draw s neighbors 
xi,...,x s of o as in the configuration model. The remaining neighbors 
y s +i, ■ ■ ■ , Um are set to be the roots of independent Galton- Watson trees with 
probability distribution p, which in particular means that Prob(A;(y/) = k) = 
Pk, I = s + 1, ... ,m. This procedure is continued ad infinitum. For s > 1, 
the obtained tree contains a subtree, comprised by an infinite number of 
distinguished rays {o, X\, x^, ■ ■ • , x n , ■ ■ ■ } with the mentioned path property 
"the neighbor of x n -\ has degree k with probability p^". The set of vertices 
of this subtree is contained in that of a Cayley tree with root o and branching 
number s. For s = l, our graph is a size-biased Galton- Watson tree, see 
pp. 407-412] or [16J, in which there is only one distinguished ray. 

In the remaining part of the article we use the following nomenclature. 
By GW(s,j)), s G No, we denote the model constructed in this subsec- 
tion. Then GW(l,p) and GW(0,p) are is the size-biased and the ordinary 
Galton- Watson trees, respectively. 



3 The thermodynamic states of the Ising model 

3.1 The thermodynamic states 

Let G = (V, E) be a general tree with root o. Given x G V, the distance 
p(o,x) is the length of the path $(o, x). For n G No, by S n we denote the 
collection of vertices x such that p(o, x) = n. Put 

n 
m=0 

and let E n be the set of all edges with both ends in V n . Let also be the 
set of {x, y} such that x G S n and y G S n+ i. 

By a n we denote the configuration of spins in V n , that is, a n = {a x : x G 
V n }. In the sequel, we suppose that 

J := sup \J X y\ < oo. (3-1) 

{x,y}£E 

The Gibbs probability distribution of configurations a n at temperature T 
and h = is the following probability measure 

Kn(vn\0 = z \ £ , exp I ^ (K X y(j x (Ty + K) + ^ K xy o x i y 

(3.2) 
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where 

K xy :=pj xy , K:=(3J, p-.= l/k B T. (3.3) 

Note that K xy can have arbitrary signs, whereas K is positive. In (|3,2p . the 
first (resp. second) summand in exp(- • • ) corresponds to the interaction of 
the spins in V n with each other (resp. with the spins £ y fixed outside V n ). 
One observes that the latter interaction involves only £ y with y G S n +i, 
which constitute the outer boundary of V n . For technical reasons, we add 
positive constants to the spin-spin interaction along each edge in E n , c.f. 
The partition function is then 

-M0 = 5Z eXP (K xy O x O y + K)+ ^ K xyVx£y\, (3.4) 

where the summation is taken over all c x — il, x G . 

Given zeV, let n z G No be such that z G S„ z , i.e., n z = p(o,z). For 
these z, n z , and for n > n z , 

(T z =±l 

(3.5) 

is the magnetization at z G V n in the state 7r n (-|£) ([3.2|) . Here, for a z = a, 
Qn,z( a \0 = ^2 n n(<7n\0, a = ( 3 - 6 ) 

T n - cr z =a 

where the summation is taken over all a n with the fixed a z = a. According 
to the theory of Gibbs states [T], the global thermodynamic state is unique 
if, for all z, the limits 

g z (a) := lim £„. 2 (a|£) (3.7) 

exist and are independent of £. In this case, by (|3.5p the global magnetization 
M z exists and 



M z = lim M HjZ (0 = Qz (l) - Q Z (-1). 

As h = 0, we have 

M n>z (0 = -M n:Z (-0, (3.8) 

and hence M z = —M z . Therefore, the uniqueness of thermodynamic states 
occurs if and only if, for all z G V and £, 

M z = lim M nz (0 = 0. (3.9) 

In this case, we say that the model is in a paramagnetic state. 

If the underlying tree G is random, then the measure (|3.2[) is also ran- 
dom, and hence (|3.9|) is a random event obeying the zero-one law. If the 
corresponding probability is one, we say that the model (jl.ip on G is in a 
paramagnetic state with probability one. 
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3.2 The main statement 

Let p and p be as in f)2.1|) and (I2.5p . For a £ (0, 1), we set 

b a = ^2 {k - s) a p k , b = Y^kp k \nk, (3.10) 

k>s+l k>2 

and also 

1 a ~ 1 s 7 4- 1 

K c =-ln -, K c =-ln^± 7 = l/a. (3.11) 

In the statement below we describe the Ising model (jl.ip on GW(s,p) 
graphs with p obeying f)2. 1 1) . 

Theorem 3.1 For s = 0, the Ising model is in a paramagnetic state when- 
ever K < K c . For s = 1, the same holds under the condition b < oo. For 
s > 2, the Ising model is in a paramagnetic state whenever b a < oo for some 
a E (0, 1) and K < min{K c ; K c }. 

In Section [5l we discuss the above statement in detail. Its proof is based on 
two lemmas proven in the next section. In the first lemma, we describe the 
model (jl.ip on a general tree, for which the corresponding quantities were 
introduced in Subsection 13.11 

Lemma 3.2 Let K > be such that 



q(K) := exp(4K) - 1< 1. (3.12) 

Then for each z £ V ; any n > n z , and arbitrary £ and rj, we have that 

\M n>z (Z) -M n ,M\ < 2[q(K)] n ^\S n \. (3.13) 

In the next lemma, we describe the model GW(s,j>) introduced in Subsec- 
tion! 



Lemma 3.3 Let L n be as in (2.1), s > 2, and b a < oo for some a £ (0, 1). 
Then, for every c obeying c > s 7 and c > a — 1, with probability one 

c~ n L n — > W £ [0, oo), as n — > oo. 



Proof of Theorem 13.11 For s = , by construction |S n | = L n is the number 
of individuals in generation n, see (j2.3|) . and a — 1 is the mean number of 
offsprings in the corresponding Galton- Watson process. Then the limit 

lim . L " (3.14) 

rw+oo (a - l) n y J 
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with probability one exists and is finite, see ()2,2p and Proposition 1.3 in [6, 
page 20]. If K < K c , then q(K) < q(K c ) = l/(a — 1) < 1, and hence we can 
apply (|3.13p . which yields that with probability one 



[g(K)(a-l)] n (a _" i)n -X), n->+oo, (3.15) 

which in turn by ()3.13|) yields (j3.9j) . and hence the proof for this case. 

For 5 = 1 , by (|2.10|) we have L n = 1, and hence |S n | = L n = 1 + L n . 
Since b < oo, with probability one we have 

£(a - l)~ n Y n < oo, (3.16) 

n>l 

see Proposition 6.2 in [61 page 50]. Then by Theorem 6.1 in [6, page 50], 
L n /(a — l) n with probability one tends to a finite limit. Then the proof 
follows as in the case of s = 0. 

For s > 2, the random variables Y n with different n are no more identi- 
cally distributed, and hence b < oo is not enough to get (|3.16p . Instead we 
use a more restrictive condition b a < oo, under which we get, see Lemma 
13.31 that c~ n |S n | = c~ n L n = c~ n L n + c~ n L n with probability one tends to 
a finite limit since c > s, cf. (|2.10p . Note that 

K(c) := i ln ^ < min{^ c ; K c }. 

Hence, for K < K(c), we have q(K) < 1/c, and the proof follows as in 
431511 . 

4 The proof of the lemmas 
4.1 The proof of Lemma 13.21 

Let G = (V, E) be a general tree. For a path by we denote the set of 
edges {xi-±,xi} of $. Since G is a tree, there exists exactly one path 
for any x and y. Recall that E n stands for the set of edges with both ends in 
V„. Thus, G ra = (V n , E n ) is a finite graph; G' = (V, E') is called a subgraph 
of G n if V C V n and E' C E n . A path & is said to be in G' if E^cE'. Two 
vertices z, x G V are said to be disconnected in G' if there is no path i?(z, x) 
in G'. 

By (|3.5|) and (|3.2|) . employing replica spins a we obtain the following 
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Meyer-like expansion 

Af n , z (0 - M n>z (ri) = ^ (a z - a z ) TT n (a n \OM^n\v) (4-1) 



kiEE (o-.-&,)r(E')* n (e > i7), 



where 



Z n {£)Z n (rj) f 

C Ltji cr n ,(7. 



r(E')= J] T *y> 

{x,y}eE' 

T xy = exp [{K xy a x a y + K) + (K xy a x a y + if)] - 1, (4.2) 

*n(£, 7?) = ] | exp [ifj-j, ((T^j + a x T]y)] . 

For simplicity, we do not indicate the dependence of the latter functions on 
the spins. One observes that each T xy > due to our choice of the spin- 
spin interactions in (|3.2p . cf. (|3.1|) . and that ^ n (£,v) depends only on a x 
and a x with x G S n . Fix some E' in the last line in (|4.ip and consider the 
subgraph G' with the edge set E' and the vertex set V n . If in G' the vertex 
z is disconnected from each x G S n , then the sums over the spins <7 z ,a z 
and over a x ,a x with i £ S„ get independent and hence the left-hand side 
of (|4.ip vanishes as the term (<r z — a z ) is antisymmetric with respect to the 
interchange a O cr, whereas all r xy are symmetric and the only break of this 
symmetry is related with the fixed boundary spins £(y) and r/(y), y G S n +i. 
Therefore, each non-vanishing term in (|4.ip corresponds to a path $(2, x) 
connecting z to some x G S n . Let us take this into account and rewrite 
(|4.ip as the sum over the subsets of E n containing the edges of at least one 
such path. Let O n (z) be the set of all paths connecting z to S n , and for 
$ G n (z), let £$ be the family of subsets of E n each of which contains E#. 
That is, S# = {E' C E n : E# C E'}. Note that, for distinct •&,■&' G & n {z), 
the corresponding families and £#/ are not disjoint - they include those 
E' which contain both E$ and E#/. Finally, set 

f= U ^> ( 4 - 3 ) 

that is, £ contains all sets of edges E' C E n such that the corresponding 
graph G' contains at least one path connecting z to some x G S n . Then 
(|4.ip takes the form 

M n , z (0 - M n , z { V ) = ^ ^ EE( ff ^ **) r ( E ')*n(e, r/). (4.4) 
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In view of the positivity of all T xy > 0, this yields 

1^,(0-^)1 < £ £ r(E0*„(M (4.5) 

Note that the sum 'Yl'de0„(z) Y^E'e£# contains the same summands as ^E'e£' 
see (|4.3|> . but a part of them are repeated, which yields the second < in (|4.5|) . 
From (|4.2p and f)3. 12|) we see that < T xy < q(K), for any x, y E V. Hence, 

r(E„) < 

which yields in (j4.5|) 



x E II ( 1+r ^) E r ( E ")^n(e,r?) 

ff»,fft.{i,!/}eE fl E"cE n \E^ 

= 2 ^ [<7(i0]U*H. 

Here we have taken into account that, see (|4.ip . 

E II ( 1 + r ^) E n^nitri) 

ff n,5n{i,!/}eE # E"CE n \E tf 

= E II {1 + r xy )^n(^r 1 )=Z n (t)Z n {r ] ). 

<? n ,d- n {x,y}&E n 

Note that a similar construction was used in deriving (2.16) and (2.17) in 
[17] . see also Lemma 3.1 in [T3]. Since n z is defined by the condition z E S Hz , 
the shortest path in n {z) has length n — n z . Furthermore, for every x E S n , 
there exists exactly one path $ E n (z). Then, for (/(if) < 1, the estimate 
fl3~T3l) follows from (fPj) . 

4.2 The proof of Lemma [3731 

If we prove that, for c > s 7 , with probability one 

5>-"Y„ <oo, (4.7) 

n>l 
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see (|2.8p . then the proof that L n /c n — > W G [0, +00), with probability one, 
follows as in the proof of Theorem 6.1 in [61 page 50]. The first summand 
of L n in (|2.7p is under control in view of the assumed property c > a — 1. 

As a S (0, 1), by the standard Minkowski inequality, cf. \18\ Theorem 8, 
page 319], as well as by (I2.10p and (|3.10p . we have 

s n — 1 

(Y:)<Y,((x 3 - s y) = s n - l b a . 

Then, for any q > c 7 , by Markov's inequality, cf. [181 page 311], 

P n := Prob {q- n Y n > l) = Prob (Y« > q an ) < ^ (^\ " . (4.8) 
On the other hand, for c > s 7 , we take g 6 (s 7 ,c) and rewrite 

E c "^ = EO n (^)' ( 4 - 9 ) 

n>l n>l 

For such g, by (|4.8p we have that Xm>i < 00 • Thus, the Borel-Cantelli 
lemma, cf . [181 P a g e 320] , yields that with probability one only finitely many 
of the events q~ n Y n > 1 may occur, which by (|4.9p readily yields (|4.7p . 



5 Concluding comments and remarks 
5.1 The graph model 

Let us look at the structure of generation n in GW(s,p), cf. ()2.7p - (|2.10p . 
For s = 1, we deal with a size-biased Galton- Watson tree, cf. [7] and |16j . 
and b < 00 is known as the X log X condition of the Kesten-Stigum theorem, 
see also [6l page 23]. If it holds, the influence of the distinguished ray on the 
structure of S n is asymptotically negligible, see Section 3 in [TB]. For s > 2, 
the number of distinguished vertices L n in generation n increases in such a 
way that the sequence {a~ n L n } n ^ with probability one tends to a random 
variable W such that Prob(W > 0) = 1 and (W) = 1. Here 

a = s - (s - l)p 2 - (s - 2)p 3 p s . 

Note that a E (1, s], and a = s if and only if P2 = ■ ■ ■ = p s = 0, and hence 
the least value of X is s. In this case, a — 1 > s. If a — 1 > s 7 , L n asymptot- 
ically 'behaves' like in the ordinary Galton- Watson tree, which means that 
the production of offsprings with probability p 'dominates' the immigration 
from the distinguished part of the population, which is always the case for 
GW(l,p), cf. (|2.2p . However, for a — 1 < s 7 , the main contribution into 
L n comes from the immigration. If o — 1 is close to one and s 7 S> 1, even 
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for big n the structure of (V„, E n ) is quite close to that of the corresponding 
subgraph of GW(0,p), i.e., as in the case of the configuration model. Note 
that s 1 is big whenever (a) s is big, and hence the distinguished subgraph is 
'big'; (b) a is small, and hence the offspring production in the distinguished 
subgraph is very intensive. 

Noteworthy, the procedure of formation of edges in GW(s,j)) resembles 
that of [HIE]. The 'value' specified in [S] by a positive ui, in our case appears 
as a mark k, distributed according to p with pk 'decreasing at infinity', cf. 
"inverse mass-action principle" of [9j. Distinguished individuals have the 
right to choose neighbors. Each chooses according to its own mark at most 
s distinguished neighbors, preferring those with big marks. The restriction 
of the number of distinguished neighbors to s corresponds to the restriction 
of the number of outgoing edges in the Cameo principle [9]. 

One observes that GW(s,j)) is a natural generalization of the size-biased 
Galton- Watson tree GW(l,j>). In a separate work, we plan to study this 
model in more detail. Another generalization of GW(l,p) can be the model 
in which the distinguished part of X is not X — {X — s)+ like in GW(s,p), 
cf. (|2.7p . but some more general increasing function of X, or is random. 
Further generalizations of this kind can be obtained in the approach of |19j . 



5.2 The thermodynamic states 

In 1989, R. Lyons [15] proved that, for any tree, the critical temperature of 
the model (jl.ip with h = and J xy = J > is 

T c = — i-, (5.1) 

ks coth g 

where g is the branching number of the underlying tree. For an infinite (al- 
most surely non-extinct) Galton- Watson tree, with probability one g equals 
the mean number of offsprings, see Proposition 6.4 in |20j. Hence, for 
GW(0,p), we have g = (k 2 ) / (k) — 1, which yields in ()5.ip 

The very same expression for T c was 'rediscovered' in 2002 in \12\ 113] . It 
is significant that Lyon's paper [15] was quoted^] in [3] and even the same 
formula (|5.ip was discussed, see the text between equations (89) and (90) 
on page 1304 in [3] or the text between equations (89) and (90) on page 33 
in larXiv:0705.0010l 

Now let us turn to the discussion of Theorem 13. 1\ which describes the 
model (jl.ip with arbitrary signal of the intensities J xy . The only condition 



5 In the preprint arXiv:0705.0010 version of [3], both papers 15, 20] were quoted. 
6 The antiferromagnetic case J xy = —J, J > 0, can be reduced to the ferromagnetic 
one by changing signs of a x , x G S2n- 
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is that they are uniformly bounded, cf. (|3.ip . As follows from the estimate 
(|3,13p . the model (jl.ip on a general tree can be in a paramagnetic state if 
the number of vertices in the sphere S n admits the control |S n | < c n for 
some c > 1. Note that in the example of [13], |S n | ~ n\. For graphs other 
than tree^j], one has to control O n (z), cf. (|4.6p . The results of Theorem 
13.11 can naturally be extended to models with continuous spins and bounded 
interaction, as e.g. in [H]. In a separate work, we plan to study the case 
of unbounded interaction intensities, which includes also random J xy , as 
e.g. in the Edwards- Anderson model. For a ferromagnetic Ising model on a 
GW(0,p), the exact value of the critical temperature is given in f)5. 1 j) with 
g = a — 1. The corresponding result of Theorem 13. II is that the model is in 
a paramagnetic state for 



It is naturally less precise as we cover the case of arbitrary signs of J xy . 

Regarding the model GW(0,p) studied in [121 [13], by Theorem 11.11 we 
get that the Ising model with arbitrary signs of J xy obeying (|3.ip is in a 
paramagnetic state if 



It is believed, cf. [22], that the Internet can well be modeled as a scale-free 
graph with A = 5/2, which corresponds to the choice of a < 1/2, see (|3.10p . 
As a byproduct of the result of Theorem 13.11 we obtain the following, cf. 



Proposition 5.1 Let 6 (0, 1) be the Bernoulli bond percolation probability 
on GW(s,p) graphs with s > 2 andp obeying A2.1\) and b a < oo. Then with 
probability one there is no giant component if 



where q(K) is given by \3.12}) . and K c and K c are as in Theorem \l.l\ 
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7 For more on phase transitions on graphs see |21| . 




a- 1 



a 




[221 [23] 



6 <mm{q(K c );q(K c )} 
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